Decoherence and localization in tunneling process under influence of one external 

degree of freedom 



O 
O 

(N 

C 



> 
in 

(N 

o 
o 

O 

^— > 

c3 

CT 1 - 



X 



Quanlin Jie, 1 '^ Bambi Hu, 2,3 and Guangjong Dong 2 

'Department of Physics, Wuhan University, Wuhan 430072, P. R. China 
2 Department of Physics and Centre for Nonlinear Studies, 
Hong Kong Baptist University, Hong Kong, China 
J Department of Physics, University of Houston, Houston TX 77204-5506 
(Dated: February 1, 2008) 

We investigate numerically the tunneling effect under influence of another particle in a double well 
system. Such influence from only one degree of freedom makes decoherence and quantum-classical 
transition, i.e., suppression of the tunneling effect. The decoherence happens even for cases that the 
influence is from a particle of very small mass, and it has virtually no effect in the corresponding 
classical dynamics. There are cases similar to dynamical localization that the suppressed tunneling 
rate is several times smaller than the classical counterpart. This result is relevant for understanding 
quantitatively the dynamical process of decoherence and quantum to classical transition. 
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Decoherence is a time honored problem that attracts 
many investigations for the mechanism of quantum to 
classical transition 1, 2J. There are growing interests in 
decoherence in the field of quantum information |J, 1^, |5j . 
A common approach to decoherence is to study sys- 
tems that thermally interact with large degrees of free- 
dom, such as enviroment. This approach needs statis- 
tical method to handle the large degrees of freedom of 
the enviroment, and one usually resorts to master equa- 
tion |y,l3)|8| or semiclassical approach [9j for the reduced 
density matrix. There is little investigation of decoher- 
ence that involves few degrees of freedom. One usually 
expects quantum effect prominent when there are only 
a few degrees of freedom involved in a physical process. 
However, as shown in Ref. |lfj| . if a time dependent sys- 
tem is classical chaotic, such as kicked rotator, decoher- 
ence happens under influence of its one internal degree of 
freedom. The chaotic nature of the underlying dynam- 
ics makes the internal degree of freedom to behave as a 
noise HI 03. 

In this Letter, we investigate numerically the tunnel- 
ing process of a particle in a one dimensional double well 
system under influence of another particle. Tunneling 
is a basic quantum effect that has far reaching applica- 
tions. Understanding decoherence in tunneling effect is 
of fundamental importance 0, 0, ^| . Our numerical 
results show that the tunneling can be totally suppressed 
under influence of only one external degree of freedom. 
Even such influence has virtually no effect in the corre- 
sponding classical dynamics, it is still able to cause deco- 
herence and hence makes the tunneling rate to approach 
the classical counterpart. Another result is that, similar 
to dynamical localization, the suppressed tunneling rates 
are several times smaller than the classical counter part 
in some cases. 

The Hamiltonian of the two-particle system reads H = 
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Vi{xi) + v(xi — X2). Here we consider 



the motion of particle 1 under influence of particle 2 
through the interaction v(x\ — £2). Our tests indicate 
that the results are rather insensitive to the exact form 
of the interaction provided that the total system is not 
integrable classically. We use an attractive interaction 
v{x\ — X2) — 57(^1 — X2 — lo) 2 in the following discus- 
sions. The double well potentials for both particles are a 
Gaussian shaped barrier at the bottom of a bounded har- 
monic potential, Vi(xi) — -^kix\ + \i exp(— x 2 /a 2 ). This 
form of double well performs similar to the common used 
form X(x 2 — a 2 ) 2 . Yet it is more flexible to control the 
height and width of the barrier by adjusting parameters 
Aj and a$, and hence is more relevant for experimental 
realization by imposing a barrier at the minimum of a 
parabolic trap jlq . 

We consider each of the two particles initially as 
a Gaussian wave packet centered at the bottom of 
one well with vanished central momentum, (f>o(xi) = 
Aexp(— a(xi — xoi) 2 /2h), where A is a normalization 
constant, a is the width parameter, and Xoi is center of 
the wave packet. When the interaction v(x\ — X2) van- 
ishes, the behavior of the particle 1 is well known [l7j |. 
The quantum motion has virtually only finite number 
of frequencies, i.e., the initial wave packet is effectively 
spanned by only finite number of eigenfunctions of the 
Hamiltonian. In our model, only 2 frequencies actually 
dominate the one dimensional tunneling processes. 

When the particle is under influence of other degrees of 
freedom, the tunneling behavior may changes completely. 
From our numerical tests, the suppression of the tunnel 
effect happens under influence of only one degree of free- 
dom, and such influence is classically negligibly small. In 
Fig. we show the tunnel effect of particle 1 under in- 
fluence of particle 2 with a mass that is I0 4 times smaller 
than that of particle 1, = 10 _4 mi. The top of Fig. 
is the tunneling rates of particle 1 versus time for various 
interaction strength 7. The particle I with unit mass is 
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FIG. 1: (Top) Tunneling rates versus time for interaction 
strength 7 = 0.01 (dot-dashed), 0.05 (dashed), and 0.2 (thick 
solid), respectively. The thin solid line is the classical coun- 
terparts. (Bottom) The von Neumann entropy versus time. 
From bottom to top, 7 = 0.01, 0.02, 0.05, 0.1, 0.2. 



initially rest on the bottom of right well. The tunneling 
rate (also known as tunneling probability) is the prob- 
ability to find the particle on the other side of the bar- 
rier, T r = J^p^x^dxi. Here pi = J \ip(xi, x 2 )\ 2 dx 2 is 
the reduced density for the particle 1 by tracing out the 
other degree of freedom from wave function ip{x\,X2)- 
The double well potential parameters (in arbitrary unit) 
are h = 0.5, Ai = 3.0, oti = 1.0, and l = 0.5 (The 
Planck constant h is set to 1). The potential for the 
particle 2 is the same as that of particle 1. We set par- 
ticle 2 to be the same Gaussian wave packet as parti- 
cle 1 that initially locates in the bottom of another well 
with vanished central momentum. The width parameter 
for both Gaussian wave packets is a = 3.0. In classical 
limit, the particle 2 with such a small mass has virtu- 
ally no effect on the motion of particle 1. The quantum 
case, however, is totally a different story. The interac- 
tion with particle 2 can completely destroy coherence of 
particle 1. As shown in Fig. ^ for very weak interaction, 
7 = 0.01, the motion of particle 1 is almost unaffected. 
It is still able to tunnel through the barrier back an forth 
almost completely. When the interaction becomes a lit- 
tle bit stronger, 7 = 0.05, the decoherence is evident. 
The maximum tunneling rate is at most 50%, and such 
maximum decreases gradually in longer time. Further 
increase of the interaction strength to 7 = 0.2, as shown 
by the thick solid line, the tunneling effect is virtually 
suppressed by the decoherence. The resultant tunneling 
rate approaches to the classical result. As a comparison, 
we show the corresponding classical counterpart by thin 
solid line. It is obtained by evolving the Liouville equa- 
tion from a classical density distribution that is identical 



to the quantum Wigner function initially. Similarly, the 
classical counterpart of the tunneling rate is the proba- 
bility to find the particle on the other side of the barrier. 

Since the mass of particle 2 is smaller than particle 
1 by several orders, its effect on particle 1 is just like 
a small weightless dust on a heavy body in the classi- 
cal limit. The classical motion of particle 1 is indeed 
quite regular, and is almost the same as the unperturbed 
one dimensional case. The classical motion of particle 
2, however, is irregular. The interaction with particle 
1 is a major force for the chaotic motion of particle 2 
which bounces drastically between the two wells. Such 
irregularity of particle 2 makes the total wave function 
of the two-particle system to be irregular. Entanglement 
between the two particles leads to the decoherence of 
the quantum motion of particle 1 . The irregular motion 
of particle 2 poses some difficulties in numerical calcu- 
lations. The wave function relates to particle 2 spreads 
widely in both coordinate and momentum space. One 
needs dense grid that covers a wide range. Our actual 
calculations employ periodic boundary condition to par- 
ticle 2. This indeed restricts particle 2 to a ring. From 
our test, this restriction does not alter the behavior of 
particle 1, but saves computation efforts remarkably. 

There are more frequencies involved in the two-particle 
system, and the structure of the frequencies is much com- 
plicated than the single particle case. One can obtain the 
spectrum density by Fourier transformation of the auto- 
correlation function. Since the mass of particle 2 is much 
smaller than that of particle 1, the energy spectrum dis- 
tributes in a wide range with band like structure. Each 
band consists of several levels. It is virtually impossible 
to have constructive interference between those frequen- 
cies on the other side of the barrier. 

The complex structure of the involved frequencies and 
the corresponding eigenfunctions relates to the entangle- 
ment of between the two particles. Such entanglement 
makes expansion of initial wave packet in terms of the 
eigenfunction to be complicated. An initially wave packet 
must evolve into entangled form in such system. The bot- 
tom of Fig. ^shows the von Neumann entropy of particle 
1 versus time for various interaction strengths. The von 
Neumann entropy, S — Tr(p\ ln(pi)), is a common used 
measure for the entanglement between two particles. By 
comparison with top part of Fig. ^for the correspondent 
interaction strength, it is evident that the decoherence 
in the tunneling process corresponds to the entanglement 
between the two particles. Increase the interaction be- 
tween the two particles leads to stronger entanglement 
(larger entropy), and hence deeper suppression of the 
quantum effects. An initially Gaussian wave packet needs 
some time to develop into the final entangled form. The 
speed of such development increases with the interaction 
strength. 

When the interaction is too strong, the two particles 
can bind together like a single one. When this hap- 
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FIG. 2: The tunneling rate versus interaction strength 7. 
The circle (diamond), down and up triangle symbols corre- 
spond to quantum (classical) mean value, typical maximum 
and minimum value of the quantum tunneling rates, respec- 
tively. (Inset) The density distribution p of particle 1 versus 
the coordinate X at time t = 400 for 7 = 0.2. 



pen, the quantum coherence emerges again. In other 
words, when the interaction strength is stronger than 
some value, the degree of freedom for the relative mo- 
tion is hard to excite, and there is enhancement of the 
quantum coherence with the increase of the parameter 7. 
Note that when two particles bind together, they behave 
like a particle of mass m\ + mi in a potential v\ + V2- 
Here mi, m,2 and v\, V2 are the masses and potentials 
of particles 1 and 2, respectively. Fig. [21 shows averaged 
tunneling rate of particle 1 versus the interaction param- 
eter 7 for both quantum and the corresponding classical 
cases. The open circles and diamonds are quantum and 
classical mean tunneling rate, respectively. The down 
and up triangles are a pair of typical maximum and min- 
imum tunneling rates chosen from a tunneling cycle. The 
range between such a pair of the triangles relates to the 
degree of coherence in the tunneling process. The initial 
states of the two particles and potential parameters are 
the same as that of Fig. 2] In classical limit, because of 
small mass, the particle 2 possesses little energy and mo- 
mentum to affect the motion of particle 1. In contrast, 
the quantum motion of particle 1 exhibits decoherence 
by interaction with particle 2 for a wide range of inter- 
action strength. As the interaction strength parameter 7 
grows from 0, the averaged quantum tunneling rate ap- 
proaches to classical limit exponentially. Then quantum 
and classical mean tunneling rates coincide for a wide 
range of strength parameter 7. When 7 > 0.8, the bind- 
ing energy between the two particles becomes hard to ex- 
cite, quantum tunneling rate increases and deviates from 
the classical counterpart. Further increase of 7 makes 
the quantum motion to gain more coherence, and to be- 
have more like the motion of a single particle. In all 
the cases, corresponding to the regular classical motion, 
the reduced density distribution for particle 1, as shown 
in the inset of Fig. [21 is quite regular. Another point 



to note is that when the quantum mean tunneling rate 
approaches to the classical counterpart, its fluctuation, 
as shown by the range between up and down triangles, 
is larger than the classical case. This is resultant from 
the fact that the dynamics of the two-particle system has 
only finite number of frequencies. 
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FIG. 3: Tunneling rate versus time for interaction strength 
7 = 0.01 (Top), and 7 = 0.1 (Bottom), respectively. The 
solid and dashed lines are quantum and classical results, re- 
spectively. 



In some cases, the fact of finite frequencies makes the 
quantum tunneling rate several times smaller than the 
classical counterpart. Such cases occur at some interac- 
tion strengths and the masses of the two particles close 
to each other. We show such a case in Fig[3 The param- 
eters for the on site potentials and interaction, as well as 
the initial states of the two particles, are the same as that 
of Fig. ^except A2 = 15. The masses of the two particles 
are same, mi — m2 = 1. As shown in the top part of Fig. 
[3J when the interaction is weak (7 = 0.01), the quantum 
coherence is almost untouched. As the strength of inter- 
action increase, quantum and classical behaviors change 
in opposite way. Quantum tunneling rate is quickly sup- 
pressed by decoherence. The classical counterpart, how- 
ever, increases with the interaction strength. This is be- 
cause the two particles are initially located on the dif- 
ferent sides of the barrier, each particle has a pull force 
on the other one. As the parameter 7 increase, they pull 
each other with larger force. This leads to increase on 
the classical motion's energy to climb over the barrier. 
When the parameter 7 = 0.1, as shown in the bottom 
of Fig. [3J the quantum tunneling rate is about 3 times 
smaller than its classical counterpart. 

This effect, similar to the dynamical localization [l8j |. 
results from destructive interferences between involved 
frequencies. It relates to two factors: decoherence from 
interaction with particle 2 and only finite frequencies in- 
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volved in the quantum motion. In classical limit, the at- 
tractive interaction between the two particles makes each 
particle gaining energy from other one for climbing over 
the barrier. The energy transformation between the two 
particles is prominent when the frequency of the relative 
motion matches the overall motion, i.e., resonance be- 
tween relative motion and overall motion. For the given 
parameters, the system is close to the classical resonance 
when 7 = 0.1. For quantum motion, however, both rel- 
ative and overall motions contain only finite frequencies. 
There is little chance that the two kinds of frequencies 
match. In other words, the quantum motion doesn't feel 
the resonance. On the other hand, the interaction with 
particle 2 destroys coherence between these involved fre- 
quencies. It is virtually impossible to make constructive 
interference on the other side of the barrier. Thus the 
quantum motion behaves like a classical motion without 
the resonance effect. This is in fact a kind of quantum 
dynamical localization to the classical resonance. 



In conclusion, we numerically show that decoherence 
and quantum-classical transition occur in the tunneling 
process under influence of only one external degree of 
freedom. Such influence can be negligibly small in the 
classical limit. Here, the correspondent classical motion 
of the one dimensional tunneling is regular when there 
is no interaction with external degree of freedom, and 
such regularity of classical motion is almost unchanged 
by such negligible interaction. However, the interaction 
makes the motion of external degree of freedom to be 
chaotic. It is such irregularity of the external degree of 
freedom that causes decoherence in the quantum motion 
by entanglement with the external degree of freedom. 
Since there is only finite number of frequencies involved 
in quantum motion, the quantum motion is unable to 
feel classical resonance. This makes quantum tunneling 
to behave similar to the dynamical localization in the 
kicked rotator, i.e., the tunneling rate is several times 
smaller than the classical counterpart in some cases. 
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FIG. 4: Tunneling rate versus the barrier height A of the par- 
ticle 2. The circle (diamond), down and up triangle symbols 
have the same meaning as that of Fig. [5] (Inset) Density 
spectrum pE versus energy E for A = 15. 

To verify the above explanation, we rerun the calcula- 
tion with different barrier hight for particle 2. The result 
is shown in Fig. 0] The symbols and their meanings are 
the same as that of Fig. except the horizontal axis is 
the barrier's height of particle 2. As the height of the bar- 
rier increase, the particle 2 becomes harder and harder 
to move across barrier, and it is virtually trapped in the 
left well. This makes the classical motion alway from 
the resonance region. When this happens, as shown in 
Fig. 01 the classical tunneling rate drops. The restriction 
on the particle 2 has virtually no effect on the quantum 
motion of particle 1. Thus when the classical dynamics 
leaves the resonance region, the quantum tunnel rate ap- 
proaches to its classical counterpart. As shown by the 
inset, when the masses of the two particles close to each 
other, there are less involved frequencies than the cases 
shown in Fig. and Fig. This means the involved 
frequencies have larger level space, i.e., the frequencies 
are more discrete. 
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